Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 1 February 2008 (MN MfeK style file vl.4) 



The role of extended corona in formation of emission lines 
and continuum in AGN 

A. Kurpiewski 1 , J. Kuraszkiewicz 2 , and B. Czerny 2 

1 Astronomical Observatory of Warsaw University, Aleje Ujazdowskie 4> 00-478 Warsaw, Poland 

2 TV. Copernicus Astronomical Centre, Bartycka 18, 00-716 Warsaw, Poland 



, 1 February 2008 

OS 



> 

O 
(N 



> 

oo 



ON 

0,; 

6 

c3 



ABSTRACT 

We study a model of an accretion disc surrounded by an extended corona of the 
temperature 10 7 — 10 s K. This corona modifies the disk spectrum since it redirects 
a significant fraction of the emission from the central parts towards the more distant 
parts of the disk. The same corona is indirectly the source of the broad emission lines 
because we expect the formation of cool clouds at the basis of the corona due to 
thermal instabilities. We model the number of the clouds and their motion through 
the corona adopting a few different physically sound assumptions. 

Comparing the predicted optical/UV continua and emission line ratios and profiles 
with the observed values we favor a particular model of a typical quasar. It radiates 
at ~ 0.5 of the Eddington luminosity, the corona surrounding the disk is additionally 
heated in excess of the Inverse Compton heating, and the broad line clouds are most 
probably destroyed when accelerated vertically above the sound speed within the 
corona. 
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1 INTRODUCTION 

There is a general agreement that the nuclear activity of 
galaxies (with exception of starburst galaxies) is powered by 
accretion of interstellar gas onto a massive black hole (e.g. 
Rees 1984). The accretion flow most probably forms a kind 
of disc structure due to an excess of angular momentum. 

The evidence of such disc-like structure is observed di- 
rectly in HST pictures of the galaxy M87 at a distance of 
several pc from the center (Harms et al. 1994) and even up 
to 0.009 pc in the case of the galaxy NGC 4258 due to the 
presence of a water maser (Greenhill et al. 1995). Closer in 
the accretion flow is clearly two-phase. The cold gas most 
probably forms a relatively flat configuration, either in a 
form of an accretion disc or blobs. This material is embeded 
in a hot medium. 

There are several observational arguments in favour of 
this general scenario. 

Disc-like geometry of the cold phase explains the lack 
of significant absorption in quasars and Seyfert 1 galaxies as 
well as the presence of the X-ray spectral features in Seyfert 
1 galaxies due to reprocessing by cool optically thick gas 
(i.e. the reflection component) covering approximately half 
of the sky for the source of X-rays (e.g. Pounds et al 1990, 
Matsuoka et al. 1990; for a review, see Mushotzky, Done & 



Pounds 1993). Additional support for the disc-like geome- 
try and Keplerian motion comes from determination of the 
shape of the K a line from ASCA data (Fabian et al. 1994). 
The reflection component is not seen in quasars (Williams 
et al. 1992) most probably due to higher ionization stage of 
the gas (e.g. Zycki et al. 1994, Zycki & Czerny 1994). 

The hot optically thin medium is required to explain the 
generation of the observed hard X-ray emission extending 
up to a few hundreds of keV (see Mushotzky et al. 1993) . It 
most probably forms a kind of corona above the cold layer of 
the gas although an alternative view was also suggested (X- 
ray emission comming from shocks formed by gas outflowing 
along the symmetry axis, e.g. Henri & Pelletier 1991). 

The heating mechanism of the corona, its structure and 
radial extension are presently unknown. 

In the innermost parts of the flow a (sometimes con- 
siderable) fraction of the total gravitational energy of ac- 
creting gas has to be dissipated in the hot corona to pro- 
vide both the X-ray bolometric luminosity and the extension 
of the emitted spectrum into high frequencies. Most proba- 
ble emission mechanism is the Compton upscattering of the 
photons emitted by the cool gas by hot thermal electrons. 
The medium is not necessarily uniform, as we observe both 
the effect of moderate Comptonization (modification of the 
high frequency tail of the big bump emission seen in soft X- 
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ray band in a number of sources, e.g. Czerny & Elvis 1987, 
Wilkes & Elvis 1987, Walter & Fink 1993) as well as sig- 
nificant (but still unsaturated) Comptonization which leads 
to formation of the hard X-ray power law. It suggests that 
hard X-ray emission is perhaps produced in hotter, maybe 
magnetically driven compact active regions (e.g. Haardt, 
Maraschi & Ghisellini 1994, Stern et al. 1995) embedded, 
or surrounded by still hot but cooler plasma. Stochastic na- 
ture of the X-ray variability supports this view (Czerny & 
Lehto 1996). 

In the outer parts of the flow the amount of energy avail- 
able is small and the radiation emitted there does not prac- 
tically contribute to the bolometric luminosity of the source 
and these regions are even less understood than the inner- 
most flow. On the other hand the formation of the Compton- 
heated corona seems to be inevitable if only the outer parts 
of the disc are not shielded from the the radiation comming 
from innermost parts (e.g. Begelman, McKee and Shields 
1983, Begelman & McKee 1983, Ostriker, McKee & Klein 
1991 - hereafter OMK, Raymond 1993). The existence of 
such a corona is actually observed in X-ray binaries (e.g. 
White & Holt 1982, McClintock et al. 1982, Fabian & Guil- 
bert 1982). The existence of such a corona in AGN may have 
very significant influence on the observed spectrum. 

The principal role of the corona surrounding the disc- 
like flow at a radius ~ 0.01 — 1 pc is not in the direct dis- 
sipation of the energy but in redirecting the radiation gen- 
erated in the inner parts towards outer parts of the disc 
by almost elastic scattering. The direct irradiation in the 
case of AGN is not efficient unless the source of radiation 
is situated high above the disc surface. The disc surface in 
AGN does not flare in its inner parts, according to the widely 
adopted description of the radiation pressure dominated disc 
by Shakura and Sunyaev (1973) and in its outer flaring parts 
does not cover more than a per cent of the sky of the central 
X-ray source (e.g. Hure et al. 1994, Siemiginowska, Czerny 
& Kostyunin 1996). 

The irradiation of the disc surface due to the scatter- 
ing by the corona has two consequences: (i) the irradiation 
modifies the optical/UV continuum emitted by the disc and 
(ii) it leads to line formation at the basis of the corona thus 
contributing significantly to the observed Broad Emission 
Lines. Both effects are the subject of our study, with the 
aim to confirm the presence of the outer corona and to con- 
strain its properties. 

The plan of this paper is following. The description of 
the model of the corona and the method of computation of 
the continuum emission of the disc and line intensities and 
profiles from clouds forming in the disk/boundary layer is 
described in Section 2. In Section 3 we present the results 
for a range of model parameters and compare them with 
observations as well as we discuss the consequences of the 
model. Conclusions are given in Section 4. 

2 THE MODEL OF THE CONTINUUM AND 
LINE EMISSION 

2.1 Structure of the corona and location of the 
source of photons 

The model of the corona above the outer parts of an accre- 
tion disc (0.01 — 1 pc) is based on the theory of the two-phase 



equilibrium studied originally by Spitzer (1978) and further 
developed in the context of AGN in a number of papers, 
starting from Krolik, McKee and Tarter (1981) (hereafter 
KMT). 

The temperature of the corona is mainly constrained 
by the inverse Compton heating and cooling by the incident 
radiation. The corona is irradiated by three different ra- 
diation sources: the high-energy central UV/X-ray source, 
radiation from the central source that has been scattered 
in corona and lower energy radiation from underlying vis- 
cous accretion disc. The central hard X-ray emission and 
the central thermal emission (for the most part UV) of the 
disk are both represented as a point like source located at 
a height Hx along the symmetry axis. Although in reality 
both emission regions are extended most of the energy is 
released within a few gravitational radii from the black hole 
which is ~ 10~ 4 or less of the the outer radius of a corona so 
the extension of the central source can be neglected. In our 
calculations we take half of the total disk luminosity as the 
value of central thermal luminosity. The proportion of total 
X-ray luminosity to central thermal luminosity is discussed 
in Sec. 2.3.1. 

Folowing KMT and OMK we denote the inverse Comp- 
ton temperature of the direct radiation (i.e. from the central 
source) by Tic- For the spectrum of quasars used by KMT 
the value of T IC is ~ 10 8 K (KMT), but allowing for more 
emission from "big blue bump" reduces Tic to ~ 10 7 K 
(Mathews & Ferland 1987, Fabian et al. 1986). The pres- 
ence of additional heating besides inverse Compton process 
does not modify the basic picture (if the heating rate is pro- 
portional to the gas density, e.g. Yaqoob 1990). However, it 
can raise the temperature of the corona above the inverse 
Compton limit, thus decoupling its value from the shape of 
the incident radiation spectrum. Therefore, we can treat the 
maximum value of the surface temperature of the corona, 
Tc, as a free parameter of the model. 

Including scattered and disc components of irradiation, 
as well as bremsstrahlung cooling process near base of the 
corona, the temperature of the corona (constant in the ver- 
tical direction) is (OMK) 

_ Tc ^dir(r) ~t~ ^scat(r) "t" F'uiscfr) (-^vzsc(r) f^C ) , 
J-cor(r) — ~p i (,-LJ 

where 

F(r) Tdir(r) ~t~ F sca t( r '} -|- F v ^ sc ^ r j (2) 

is total radiative flux at the base of the corona from all 
sources of irradiation. Indices 'dir', 'scat' and 'vise' are re- 
lated to direct, scattered and disc components respectively 
and T v i sc is the temperature of radiation emitted locally 
from the disc (F visc ). 

The fraction of the central direct radiation (Fdir) and 
radiation scattered by the corona (F sca t) towards the disc 
at a radius r is given by simple analytical formulae of OMK 
in the case of sources with low ratio of the luminosity to the 
Eddington luminosity. 

If the luminosity of the source is closer to the Edding- 
ton luminosity the corona is optically thicker and multiple 
scatterings play a significant role. In this case the numerical 
computations are necessary but the results are given in the 
paper of Murray et al. (1994) - hereafter MCKM - for a few 
sets of model parameters. 
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We use both the papers of OMK and MCKM to cal- 
culate the direct and scattered radiative flux irradiated the 
disk at a given radius. Therefore we follow the assumptions 
about the corona structure and the geometry made in those 
papers. 

The corona extends up to 0.2 of the true maximum 
radius where the thermal energy of the gas particle ~ kTc 
is balanced by gravitational energy of the particle Gh ^ 



rc = 



GMfj, 
kT c 



10 10 

Tcs 



M 



cm, 



(3) 



where Tcs is the maximum surface temperature of the 
corona expressed in units of 10 s K, /i is the mean mass per 
particle (fi = 0.61m p for fully ionized gas of cosmic abun- 
dances) and M is the mass of the black hole. 

The ionization parameter H (KMT) at the basis of the 
corona is equal (McKee & Begelman 1990) 



=, b(r) « l-3T co 3 r8 2 (r) , (4) 

assuming that at higher densities bremsstrahlung is the only 
atomic cooling process. Since the ionization parameter is 
defined as a ratio of the incident radiation pressure to the 
gas pressure, its value determines the density at the basis of 
the corona (see Section 2.3.2). 



2.2 Computation of the optical/UV continuum 

2.2.1 Heating by viscous dissipation 

Accretion disc is heated by dissipation of the gravitational 
energy of the gas through viscous forces. The radiative flux 
which corresponds with this energy for the disc with nonro- 
tating black hole is given by (Page & Thorne 1974) 



visc(r) 



3GMM 



r) 



f(r) = 



where 



2GM 



r*(r- §) 



3.01 • lQ-°T C src 



is the Schwarzschild radius, and 

Lvisc 



M 



ec z 



(5) 



(6) 



(7) 



(8) 



is the accretion rate. The efficiency e is ~ 5.6% for a non- 
rotating black hole. In equations (5) and (6) the radius r is 
expressed in r g units. 

This picture does not leave any room for the X-ray emis- 
sion. Actually, a fraction of energy due to accretion is dissi- 
pated in the form of X-rays. However, reliable predictions of 
that fraction as a function of radius are not available (see e.g. 
Witt et al. 1996). On the other hand the efficiency of accre- 
tion is most probably higher than adopted as the black hole 
may well be rotating, increasing the bolometric luminosity 
easily by a factor few. Therefore we describe the accretion 
disc emission as above but we allow that the total lumi- 
nosity of the central source L may be higher than the disc 
luminosity L visc due to the contribution from unaccounted 
X-ray emission. 



2.2.2 Direct radiative heating of the disc 

The additional radiation flux heating the disc is the flux of 
radiation from the central UV/X-ray source with luminosity 
L 

L 



- dir(r) 



(l-A) 



47r(r 2 + H%) 



(9) 



where A is the albedo of the disc, 8dir(r) is the angle between 
the incident ray and the normal to the disc at r, and 



/ = 



L/Anr 2 



(10) 



is used by OMK and MCKM dimensionless factor, which is 
the ratio of the radiation flux at the base of the corona to 
the unattenuated flux from the central source (in this defi- 
nition in the denominator the component H% is neglected). 
To determine the factor f dir we use the results of OMK 
or MCKM papers as we mentioned in Sec. 2.1. We assume 
A=0.5 as an appropriate value for quasars because the disc 
surface is partially ionized (see e.g. Zycki et al. 1994). 

2.2.3 Heating by corona 

To describe the effect of irradiation by photons from the cen- 
tral source scattered towards the disc surface by extended 
corona we also use the method described in Sec. 2.1 to de- 
termine the scattered radiation flux 

Fscat(r) = (1 — A) fscat(r){cOs0 SC at), (H) 

where {8 sca t} is the angle between the direction from the 
scattering point and the normal to the disc at r, averaged 
over the whole volume of the corona. We determine the 
factor fscat(r){cos9 aca t) from OMK or MCKM papers. We 
adopt the same albedo as for the direct flux. 



2.2.4 Computation of the spectrum 

The effective temperature of the disc photosphere is given 
by 



visc(r) 



H - -*air(r) ^scat 



(12) 



Assuming the black body emission from the disc, one can 
calculate the shape of continuum as follows 



r ma . 

U = 2n J 



rB v [T eJf(r) ]dr, 



(13) 



where B v is the Planck function and r max is the outer radius 
of the disc (in calculations we assume r max — rc). 

We neglect the modification of the disc spectrum due 
to electron scattering. Although such effects were discussed 
by a number of authors (e.g. Czerny & Elvis 1987, Ross & 
Fabian 1993, Shimura & Takahara 1995) these corrections 
strongly depend on the accuracy of the description of atomic 
processes as well as assumptions on the disc viscosity. 



2.3 The spectrum of the central source and the 
emission lines 
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2.3.1 Incident radiation spectrum 

The computations of the continuum given in Sect. 2.2 cover 
only the optical/UV band as only the knowledge of the to- 
tal X-ray luminosity, but not the shape of the primary X- 
ray emission, was required to compute the thermal emission 
of the disc dominating in this spectral band. On the other 
hand the computations of the strength of the emission lines 
require the determination of the spectrum in the EUV and 
X-ray band as well. 

Therefore, for the purpose of calculating emission lines 
we use the spectral shapes which approximate well the ob- 
served overall spectra of Seyfert galaxies and quasars. 

The shape of the big blue bump is parametrized in a 
similar way as by Mathews & Ferland (1987). However, we 
adjusted parametrization in agreement with the present ob- 
servational data and we adopted interpolation to give the 
ratio of the bolometric luminosities of the big blue bump to 
X-rays equal ~ 1 for Seyfert galaxies and ~ 10 for quasars. 
The details of this parametrization are given in Appendix. 
The inverse Compton temperatures for these two spectra are 
5.8 • 10 7 K for Seyfert galaxies and 1.1 • 10 7 K for quasars. 



2.3.2 Local line emissivity 

We assume that the high ionization emission lines come from 
the clouds which form in the narrow intermediate zone be- 
tween an accretion disc and a hot corona. Thermal insta- 
bility of irradiated gas at intermediate temperatures causes 
discontinuous transition between a disc and a corona (Begel- 
man, McKee & Shields 1983), if only radiative processes are 
taken into account. In realistic situation, when some level of 
a turbulence is present in the medium, we may expect the 
spontaneous formation of the cool clumps embedded in the 
hot coronal plasma in a relatively narrow transitory zone 
(Rozahska & Czerny 1996). 

The column density of the clouds is not determined pre- 
cisely but the estimates based of the relative efficiency of 
condution and thermal processes give values of the same or- 
der as values expected for the clouds forming the Broad Line 
Region. Therefore we assume for all the clouds the column 
density equal 10 23 cm~ 2 . 

Cloud parameters (density n c i and temperature T c {) at 
a given radius are determined by two requirements. This first 
condition is the pressure equilibrium with the hot medium 
at the basis of the corona. The second condition, in the case 
of optically very thin plasma, should be the settlement on 
the lower stable branch at the E — T curve (KMT). How- 
ever, this branch is not applicable for media of higher optical 
depth. Since we do not expect clouds to cool below the local 
effective temperature of the disk surface, we assume their 
temperature to be at that value, T ct = T e ff. Therefore we 
adopt the condition 



T, 



n c l(r) — »1-d(r) — W 6 ( r ) 



cor(r) 



T, 



ff(r) 



nb(r) 



cor(r) 
T c l( r ) 



(14) 



where n 6 ( r ) is the density at the basis of the corona deter- 
mined from the value of the ionization parameter H& and 
Tcor(r) and T e ff(r) are given by eq. (1) and (12). 

We discuss the kinematics of the clouds in Sect. 2.3.3 
since it is essential for computation of the line profiles. How- 
ever, cloud motion influences the line intensities as well. 



Since the clouds are blown out radiatively from the for- 
mation region they do not form a flat layer on the top of a 
disc but they are fully exposed to the incident radiation flux 
so the inclination angle of the direct incident flux does not 
have to be included through the cosine factor, as it is the 
case for the disk surface. 

We calculate the emissivity of several emission lines 
(see Table 5 for the list) using the photoionization code 
CLOUDY. The calculations are made for a grid of radii re- 
sulting from the used range of n c i : 10 9 — 10 13 cm -3 (the same 
for all models). For each radius separately we calculate the 
contribution to emission lines assuming the local value of 
the total heating flux F (r) given by eq. (2), (5), (6), (9), (11), 
density n d ( r ) from eq. (14) and fixing the column density 
N H = 10 23 cm -2 . 



2.3.3 Radial dependence of number of clouds 

Clouds forming in the transition layer between the disk and 
the corona does not necessarily cover the disk surface uni- 
formly. Local number of clouds give weight to the local emis- 
sivity thus influencing both the line ratios and line profiles. 
The number of clouds, Nr r \, existing at a given radius de- 
pends both on the cloud formation rate, N( r ), and expected 
life time of a cloud, i( r ) 

N (r) = N (r) t (r) . (15) 

As the detailed process of cloud formation and destruc- 
tion is not well understood we discuss a few representative 
cases based on available estimates. 

We consider two cases of the cloud formation process. 
In the first case we assume that only one cloud can form 
at a given moment and a given radius and the formation 
time is given by the characteristic isobaric cooling time, T( r ) 
(McKee & Begelman 1990) 

1 



case(I) TV, 



M 



n b ( r ). 



(16) 



In the second case we assume that the clouds form in the 
entire instability zone therefore the number of the clouds 
forming at the same time is related to ratio of the zone 
geometrical thickness AZ , which is of order of the Field 
lenght (e.g. Rozahska & Czerny 1996) to the size of the 
cloud, r c;(r ) = N H /n cl ( r ) 



case(II) N t 



AZ, 



(r) 



(r) 



■ n c l( r ). 



(17) 



r c l(r)T(r) 

We describe the destruction process using three differ- 
ent approaches. In the case (i) we assume that the clouds sur- 
vive only within the instability zone. As they move through 
the zone upwards under the influence of the radiation from 
the disk surface their life time is given by the travel time 
through the zone AZ. We assume that the radiation accel- 
eration is constant. It gives us the relation 

case(i) i( r ) ~ \ (18) 

In the case (ii) we assume that clouds survive even out- 
side the instability zone but undergo the destruction due to 
conduction, evaporating into surrounding hot corona. The 
timescale of such a process is given by McKee & Begelman 
(1990) 
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Table 1. The models of radial dependence of number of clouds. 



Model 


Cases from 
Sec. 2.3.3 


The radial 
dependence 


a 


(I)(i) 




b 


(I)(ii) 


const. 


c 


CO(iii) 


n b(r) 


d 


(H)(i) 


y/nb(r)/ T cl(r) 


e 


(H)(ii) 


Tcl(r) 


f 


(II)(iii) 


n cl{r) 



case(ii) t, r) ~ . (19) 

n b(r) 

In the case (iii) we assume that clouds are accelerated so 
efficiently by the radiation pressure that they reach soon the 
velocities exceeding the local sound velocity in the corona. 
Reaching this terminal velocity of ~ 2000 km/s they are de- 
stroyed by dynamical instabilities. Since we assume that the 
radiation acceleration of the clouds is independent from the 
disk radius the time needed to reach this terminal velocity is 
also constant if small variations of the corona temperature 
with the radius are ignored. In that case 

case(iii) i( r ) = const. (20) 

We summarize the six models of the number of clouds 7V( r ) 
in Table 1. 

2.3.4 Line profiles 

We assume that the velocity field of the emitting clouds con- 
sists of the Keplerian orbital motion and the outflow perpen- 
dicular to the disc surface. 

In our computations of the orbital motion we follow 
the method of Chen and Halpern (1989), including the rela- 
tivistic effects. Their procedure requires the specific intensity 
from the disc surface 

J (^e) = ' \ 27T y/2 a [ergss^cm^sr^Hz- 1 ] (21) 

where is the local emissivity of the clouds (see Sec. 2.3.2) 
and the exponential factor is related to local broadening of 
the emission line. The broadening (due to electron scattering 
or turbulence e.g.) is characterized by quantity o, and 
uo are the emitted and rest frequencies. After Chen and 
Halpern (1989) we use in our calculations cr/^o=0.03. 

The vertical motion is computed assuming that the ac- 
celeration of the clouds is constant, i.e. the velocity increases 
linearly with the distance from the disc surface. Maximum 
velocity is constrained by the cloud destruction. 

Since the line profiles are not strongly dependent on the 
detailed description of the vertical motion as the orbital ve- 
locity is usually much higher than the vertical one we do 
not distinguish between the three cases introduced in Sec- 
tion 2.3.3 and in all cases we assume that the clouds move 
under the influence of constant radiative acceleration and 
reach the same terminal velocity equal 2000 km/s, indepen- 
dently from the disk radius, as expected in case (iii). 



3 RESULTS AND DISCUSSION 

We calculate the local radiation flux at the disc surface for 
the following parameters 

M = 10 8 M Q 

s = 0.056 

T c = W S K ; 10 7 K 

L = 0.01L b ;0.34L e ;0.59Lb 

H x = 3.32r 9 ;33.2r 9 . 

The second value of the height could only be used in the 
calculations of the low luminosity cases as the numerical so- 
lutions of the high luminosity corona structure (MCKM) 
are available only for the first value. In agreement with 
generally accepted trend that low ratios of luminosity to 
the Eddington luminosity are appropriate for Seyfert galax- 
ies and luminosities closer to the Eddington limit are ap- 
propriate for quasars, we use Seyfert galaxies spectrum for 
L = 0.01Z/E and quasars spectrum for L = 0.34Le, 0.59Lb 
(see Appendix). 

3.1 IR/optical/UV continuum 

We calculate the accretion disk spectra taking into account 
the direct irradiation by the central source as well as the 
irradiation by the flux scattered in the corona, as described 
in Sect. 2.2. Figure 1 illustrates our results in loguf u — logv 
form. 

In the case of low L/LEdd ratio, the corona is never 
strong and the direct irradiation dominates, independently 
from the corona temperature. Therefore the adopted loca- 
tion of the irradiating source is of significant importance. 
No corona influence is seen when the temperature is lower, 
some redistribution of the flux towards the outer part of the 
disk is present when the temperature is higher. 

In the case of high L/LEdd ratio, the effect of the corona 
is essential and direct irradiation negligible which means 
that the location of the central source is not essential any 
longer (unless Hx would be very high indeed). The effect of 
irradiation is significant even for lower corona temperature 
and it is particularly strong for high temperature, leading 
to significant enhancement of spectra in IR/optical band. 

In order to compare the derived spectra and observed 
continuum shapes of accretion disc we calculate the spectral 
indices a (f v ~ v a ) in optical range (u ~ 10 14 ' 5 ~ 10 15 Hz) 
and in the following pieces of UV range v ~ 10 15 ' 12 10 15 32 
Hz (UV-1) and v ~ 10 15 ' 22 -=- 10 15 ' 36 Hz (UV-2). They are 
presented in Table 2. The mean values of spectral index for 
these three ranges obtained from observations by various 
groups for various samples of quasars we show in Table 3. It 
could be supplemented by the UV slope (between 15.13 and 
15.45) of the composite radio quiet quasar spectrum given 
by Zheng et. al (1996) equal -0.86, but their formal error 
does not reflect the dispersion of the contributing spectra. 
Unfortunately, equaly reliable data for Seyfert galaxies are 
not available since in the case of weaker active galactic nu- 
clei the determination of the spectral slope is complicated 
by strong contamination of the spectra by circumnuclear 
starlight. 

Comparing directly the predicted spectral slopes of the 
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Figure 1. The spectra of irradiated accretion disks parametrized by the L/L E dd ratio and the corona temperature. The dotted curves 
show the spectra of non-irradiated discs and the dashed ones illustrate the cases with Hx = 33.2r 9 . 



Table 2. Spectral indices in the selected ranges of continuum. 



Model 






auv-i 


O-UV-1 


L = 0.01L E ;T = 
H x = 33.2r 9 


10 7 K; 


- 0.5 


- 1.1 


- 1.22 


L = 0.01L B ;T = 
H x = 3.32r 9 


10 7 K; 


0.01 


- 0.74 


- 0.94 


L = 0.01L B ;T = 
H x = 33.2r 9 


W 8 K; 


- 0.06 


- 0.75 


- 0.94 


L = 0.01L E ;T = 
H x = 3.32r 9 


W 8 K; 


- 0.1 


- 0.9 


- 1.14 


L = O.UL E ;T = 


10 7 K 


- 0.19 


- 0.27 


- 0.32 


L = 0.34L B ;T = 


10 8 K 


- 0.55 


- 0.36 


- 0.4 


L = 0.59L B ;T = 


10 7 K 


- 0.32 


- 0.3 


- 0.34 


L = 0.59L B ;T = 


10 S K 


- 0.66 


- 0.32 


- 0.35 



Table 3. The mean values with standard deviations of spectral 
index obtained from observations. 



a op t 


ceuv-l 


Q C/V-2 


Neugebaucr ct al. 


Baldwin et al. 


Francis et al. 


1987 


1989 


1992 


-0.4=f0.4 


-0.91=f0.34 


-0.67=f0.5 



in Neugebauer et al. paper the emission in 'small blue bump ' is 
subtracted. 

mass would give flatter optical spectra and steeper UV spec- 
tra due to the the decrease of the disk temperature for given 
L/Lsdd ratio. It might therefore slightly favor the models 
with higher corona temperature. Unfortunately, the avail- 
able corona models are only for 10 s Mq (see Sect. 2.1) so we 
cannot support this conclusion quantitatively. 



high luminosity models (roughly adjusted to quasar lumi- 
nosities) with the observed values we conclude that they 
roughly correspond to the data, taking into account large 
errors. However, the agreement with mean values is far from 
perfect since the UV slopes are too flat. This is clearly caused 
by adopting the value of the black hole mass equal 10 8 M Q , 
actually too low by a factor 3 to 10. Higher values of the 



3.2 Emission lines 

We can easily see that the emission of the broad emission 
lines by clouds formed close to the basis of the hot ac- 
cretion disk corona is a reasonable assumption. The typ- 
ical widths of lines reaching FWHM^2000 - 10000 km/s 
correspond to Keplerian velocities over the radii range 

1.3 • 10 16 cm < r < 3.4 • 10 17 cm in the case of an accretion 
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T=10 7 K T=10 8 K 




4 i i I i i i i I i i i i I i i I i i i i I i i i i I i i '. i i 
2 3 4 2 3 4 5 
log r [r s ] log r [r g ] 



Figure 2. The examples of the radial distribution of the local 
line emissivity of four strongest lines as defined in Sect. 2.3.2 
for L/L^dd =0.59 and the corona temperature 10 7 K (left panel) 
and 10 8 K (right panel). The solid, long-dashed, short-dashed and 
dotted curves represent Ly a , CIV, Hell and NV lines respectively. 

disc around 10 s Mq black hole. It is the region covered by 
the corona and the cloud emissivity is large in this range 
(see Fig. 2). Also the line ratios should be reasonable as the 
cloud formation is based on special requirements as for the 
ionization parameter. 

Detailed predictions of the model can support this sce- 
nario and allow to put some constraints on the L/L E dd ratio, 
corona temperature and the radial distribution of clouds. 

3.2.1 Line profiles 

We present the calculated profiles for Civ and Lya lines 
only, because these lines are statistical investigated best of 
all. 

Line profiles determined by the model depend on the 
disk/corona model and on the adopted radial distribution of 
the number of clouds, reflecting various assumption about 
their formation and destruction. 

An example of the radial distribution of the local emis- 
sivity of several lines is shown in Fig. 2. It is further com- 
bined with the radial distribution of the number of clouds. 
In the case of model (b) (see Table 1) this distribution is 
preserved whilst in other models it includes additional local 
weight and either outer or inner parts of the distribution are 
enhanced. 

We calculate the profiles for four values of the inclina- 
tion angle of accretion disc: 0°, 30°, 60° and 80°. Accord- 
ing to the unified scheme of AGN supported by a number 
of strong observational data (e.g. Antonucci 1993) objects 
viewed at large inclination angles are obscured by molecu- 
lar/dusty torus and are not identified as quasars so the mean 
inclination angle is expected to be smaller than 60°. 

We compare the profiles with the mean values and ob- 
served ranges of the these two line widths. 

The first sample are exceptionally well studied line pro- 



files of quasars by Brotherton et al. (1994). These autors 
analysed two quasar samples (a small sample with high qual- 
ity UV spectra and the other from Large Bright Quasar Sur- 
vay; hereafter ALS and LBQS) and were able to decompose 
the contribution to broad lines from Intermediate Line Re- 
gion (distant and most probably spherical) and Very Broad 
Line Region. Here we use the ALS data. 

The second sample (Bright Quasar Sample, BQS) is 
from Baldwin et al. (1989) and the third one quasars from 
Wills et al. (1993) which are not decomposed into ILR and 
VBLR, as above. In our model we expect the contribution 
from broad range of radii and such a decomposition might 
not be necessary. 

In Table 4 we give the observational constraits and the 
results from our most promising models, i.e. case (a), (b) 
and (f). 

Radial cloud distribution (c) gives unreasonably broad 
profiles for both lines since they enhance too much the emis- 
sion from inner radii (the number of clouds decreases with 
radius far too fast). Models from families (d) and (e) gives 
exactly the opposite trend, leading to unreasonably narrow 
profiles. 

Model (b) also gives too narrow profiles, although not 
as narrow as case (c). Slight decrease of the cloud number 
with radius is clearly favored. 

Second and third among the class (a) models well rep- 
resent the CIV distribution and require the dusty torus to 
cover inclination angles above ~ 60°. However, these solu- 
tions do not satisfy the observational requirement that the 
Ly a is broader than the CIV. 

We conclude that the most favorable representation is 
given by model (f), L/LEdd = 0.59 and the corona tem- 
perature 10 8 K, with the torus shielding the view above 
~ 60 — 70° . This last quantity is very sensitive to the corona 
temperature. We find the solution with low temperature 
equally satisfactory if the opening angle of the torus is as 
small as 30°. 

3.2.2 Line intensities 

We compute line ratios for three models of the radial cloud 
distribution which were most promissing from the point of 
view of Ly a and CIV profiles varying also other model pa- 
rameters. The results are presented in Table 5. 

We compare them with the line ratios for quasars de- 
termined by Brotherton et al. (1994) for ALS and LBQS 
samples, Baldwin et al. (1989), Francis et al. (1991) and 
Zheng et al. (1996). 

The Ly a to CIV ratio predicted by models with the 
luminosity close to the Eddington luminosity and the mean 
quasar spectrum fall reasonably close to the observed values, 
although many models with low corona temperature tend 
to underproduce this ratio whilst high temperture corona 
models sometimes overproduce it. 

The dispersion between different observational results is 
considerable but we see that low corona temperature, con- 
sistent with the inverse Compton temperature of the quasar 
spectra, is allowed only in case (f) and L/L E dd ratio equal 
0.34; all other models require higher temperature to produce 
the required amount of Ly a . 

Model (f), favored by the study of the line profiles, with 
the corona temperature somewhat below 10 8 K seems quite 
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Table 4. FWHMs in km/s of Ly a and CIV lines profiles 



Model 






i = 


0° 


i = 30° 


i 


60° 


i = 


80° 










CIV 


Ly a 


CIV 




CIV 


Ly a 


CIV 


Model a 


L = 0.34Le;T = 


10 7 K 




3130 


3240 


7070 


7590 


■ 


12080 


12450 


13080 


L = 0.34L B ;T = 


W 8 K 




2840 


2760 


5890 


6700 


9020 


10800 


9950 


12150 


L = 0.59L E ;T = 


10 7 K 




2980 


2690 


5750 


5750 


8880 


9060 


9760 


10300 


L = 0.59L E ;T = 


10 8 K 




2690 


2690 


5194 


4900 


7710 


7400 


8440 


8070 


Model b 


L = 0.01Le;T = 


10 s K;H X = 


33.2r g 


2700 


2700 


3980 


4090 


5670 


5780 


6220 


6340 


L = 0.01L E ;T = 


W 8 K;H X = 


3.32r 9 


2800 


2800 


4460 


4720 


6670 


7040 


7150 


7900 


L = 0.34L S ;T = 


10 7 K 




2620 


2610 


3430 


3390 


4680 


4670 


5200 


5270 


L = 0.34L B ;T = 


W 8 K 




2620 


2650 


3510 


3790 


5450 


5670 


6010 


6340 


L = 0.59L E ;T = 


10 7 K 




2620 


2690 


3200 


3280 


4346 


4462 


5000 


5100 


L = 0.59L E ;T = 


10 8 K 




2600 


2610 


3790 


3790 


5670 


5780 


6300 


6450 


Model f 


L = 0.34L B ;T = 


10 7 K 




5120 


2910 


16690 


11860 


> 20000 


~ 20000 






L = 0.34L E ;T = 


10 8 K 




2870 


2870 


9320 


9770 


15360 


16020 


17240 


17970 


L = 0.59L E ;T = 


10 7 K 




4050 


2840 


11050 


9360 


16830 


15210 


19200 


16700 


L = 0.59L E ;T = 


W 8 K 




2730 


2730 


6810 


6260 


10610 


9500 


11790 


10550 


Observations 


Brotherton ct al. (1994)ALS 






Ly a : 8000 ± 600 






CIV: 6800 ± 300 




Baldwin ct al. (1989)BQS 








CIV: 2500 -7- 


S000; mean= 


5680 ± 390 






Wills ct al. (1993) 








CIV: 1970 + 10400; mean= 


=4870 ± 180 







attractive if compared, for example, with Zheng et al. (1996). 
It could reproduce well the Ly a to CIV ratio as well as 
the amount of NV and (perhaps) SilV. However, it overpro- 
duces Hell line, although this discrepancy is not so strong 
if other samples are considered. It also underproduce some 
other weak lines like A1III, Silll] and CTII]. 

However, it is possible that there is a contribution to 
the line emission from collisional excitation. In the case of 
clouds moving across the corona it is only natural to expect 
the formation of the shock at the cloud front and therefore 
some additional ionization. This effect was neglected in out 
study. 



of the contribution from collisionally heated sides of moving 
clouds may be essential. 
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4 CONCLUSIONS 

Presented results support the following view of quasars. 

Quasars are radiating at L/L E dd ratio ~ 0.5, accretion 
disk in these objects is surrounded by a hot corona with 
temperature higher than the inverse Compton temperature. 
Broad emission lines observed in their spectra come from 
clouds which form continuously at the basis of the corona 
due to thermal instabilities and are being blown out by the 
radiation pressure untill they are destroyed when they reach 
supersonic velocities, ~ 2000km/s. 

However, further studies are necessary to confirm this 
picture. For example, extension of the computations towards 
higher values of the mass of the black hole and the inclusion 



APPENDIX 

We parametrize the continua emitted by the central regions 
of the accretion flow by adopting fixed values of the energy 
index a in a number of energy bands given in logv. In the 
case of Sefert galaxies we assume 10.5 - 14.83 (a = 2. 5), 14. 83 

- 15.76 (-0.5), 15.76 - 16.12 (-1), 16.12 - 16.60 (-3), 16.60 - 
18.70 (-0.7), 18.70 - 19.37 (-0.9), 19.37 - 22.37 (-1.67). In 
the case of quasars we assume 10.5 - 14.83 (a=2.5), 14.83 

- 15.76 (-0.5), 15.76 - 16.12 (-1), 16.12 - 17.05 (-3), 17.05 - 
19.37 (-0.7), 19.37 - 22.37 (-1.67). 
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Table 5. Line Intensity Ratios 



Model 




Ly a 


NV 


CII 


SilV 


OIV] 


CIV 


Hell 


OIII] 


A1III 


Silll] 


CTII] 






1216 


1240 


1335 


1397 


1402 


1549 


1640 


1663 


1859 


1892 


1909 


Model a 


L = 0.01L B ;T = 10 8 K;H X = 


33.2r 9 


80 


19 





16 


8 


100 


18 


4 


2 


2 


3 


L = 0.01L B ;T = 10 8 K;H X = 


3.32r 9 


92 


23 





10 


10 


100 


15 


5 


1 


2 


5 


L = O.UL E ;T = 10 7 K 




114 


15 








4 


100 


76 














L = 0.34L B ;T = 10 8 X 




239 


33 





19 


8 


100 


45 


3 


1 


3 


3 


L = 0.59L B ;T = 10 7 K 




144 


20 





1 


1 


100 


58 


2 








1 


L = 0.59L B ;T = 10 8 K 




239 


25 





19 


7 


100 


42 


2 


1 


3 


2 


Model b 


L = 0.01L B ;T = W S K;H X = 


33.2r 9 


83 


4 





7 


3 


100 


6 


7 


1 


3 


14 


L = 0.01L B ;T = 10 8 K;H X = 


3.32r 9 


84 


6 





6 


4 


100 


6 


7 


1 


2 


14 


L = 0.34L B ;T = 10 7 K 




132 


5 





1 


2 


100 


24 


1 











L = 0.34L B ;T = 10 8 .Fs: 




276 


8 


1 


14 


5 


100 


21 


8 


1 


10 


14 


L = 0.59L B ;T = 10 7 K 




125 


6 





1 


3 


100 


12 


2 








1 


L = 0.59L B ;T = 10 8 K 




226 


7 





13 


4 


100 


20 


3 


1 


8 


5 


Model f 


L = 0.34L B ;T = lO 7 .^ 




210 


34 








3 


100 


144 














L = 0.34L B ;T = 10 8 X 




234 


46 





21 


9 


100 


55 


1 


2 


2 


1 


L = 0.59L B ;T = 10 7 X 




166 


25 





1 


1 


100 


83 


1 











L = 0.59L B ;T = 10 s .FsT 




250 


37 





22 


9 


100 


54 


2 


2 


2 


1 


Observations 


Brotherton ct al. (1994)ALS 




219 


111 


6 


42 




100 


31 





17 


<5 


41 


Brotherton ct al. (1994)LBQS 




154 


137 


6 


55 




100 


40 





22 




59 


Baldwin ct al. (1989)BQS 




236 


106 








100 




30 


16 




43 


Francis ct al. (1991)LBQS 




159 




4 


30 




100 




29 


46(A1III + 


CIII]) 


Zheng et al. (1996) 




192 


27 


1 


16 




100 






6 




23 
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